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ABSTRACT 

The study of the dominant population of high redshift IR-luminous galaxies 
(10^^ - 10^^ L0 at 1 < 2 < 3), requires observation of sources at the ~ 0.1 
mJy level in the mid-IR. We present the deepest spectra taken to date with the 
Infrared Spectrograph (IRS) on the Spitzer Space Telescope. We targeted two 
faint (/24 ~ 0.15 mJy) sources in the Southern GOODS field at z = 1.09 and 
z = 2.69. Spectra of the lower redshift target were taken in the observed-frame 
8-21 micron range, while the spectrum of the higher redshift target covered 21- 
37 microns. We also present the spectra of two secondary sources within the slit. 
We detect strong PAH emission in all four targets, and compare the spectra to 
those of local galaxies observed by the IRS. The 2; = 1.09 source appears to be a 
typical, star-formation dominated IR-luminous galaxy, while the z = 2.69 source 
is a composite source with strong star formation and a prominent AGN. The 
IRAC colors of this source show no evidence of rest-frame near-infrared stellar 
photospheric emission. We demonstrate that an AGN which contributes only 
a small (~ 10%) fraction of the bolometric luminosity can produce enough hot 
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dust emission to overwhelm the near-IR photospheric emission from stars. Such 
sources would be excluded from photometric surveys which rely on the near-IR 
bump to identify starbursts, leading to an underestimate of the star formation 
rate density. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
high-redshift — 



Introduction 



More than 70% of present d ay stars were forrn ed at redshifts > 0.5 in galaxies undergo- 
ing an infrared luminous phase (lElbaz et al.ll2002l ). identifiable as LIRGs (10^^ <Lir< 10^^ 
Lq) and ULIRGs (I/ir> 10 -^^ Lf;^). These source s are responsible for the peak of the faint 
24 micron number counts (jPapovich et al.l l2004l : iMarleau et al.l 120041 ) . LIRGs are partic 



ularly important at z < 1.5, beyond which ULIRGs may dominate (ILagache et al. 



2004 



Charv et al.l 120041 : IPerez-Gonzalez et all bood : Ihe Floc'h et~aD boOsL As much as 85% o f 



their luminosity may be generated by star formation (IBell et al.l l2005l : iBrand et al.ll2006l ). 
although that fraction is still uncertain. 

IR luminous galaxies are rare in the present epoch, even though they may represent a 
generic phase typical of star formation activity in distant galaxies. Loca lly, 95% of ULIRGs 
and ~ 50% of LIRGs are driven by mergers (e.g. iBushouse et al.l 120021 ). which provide the 
primary means of triggering the required amounts of star formation and AGN activity. At 
higher redshifts, however , more gas rich systems may produce sta rs at high rates under less 



extreme conditions (e.g. iMarcillac et al.ll2006l: ISaiina et a. 



the optical morphologies of ^ ~ 1 LIRGs jBell et al' 2005), which span a range of Hubble 



20061 ). This evolution is seen in 



types, from ellipticals to irregulars. 

High redshift IR luminous galaxies are easily detected by Spitzer Space Telescope imaging 



surveys such as the Great Observatori es Origins Deep Survey (GOODS; iGiavalisco et al. 



2004J : [Dickinson fc GOODS TeamI 120041) and the SIRTF Wide-Area Infrared Extragalactic 
Survey f SWIRE: Eonsdale et al.l 120031 ). These sources have been beyond the reach of most 
large Spitzer spectroscopy programs of high redshift sources, which have been limited to 



> 0.7 mJy and IRS exposures of ^ 1 hour (IHouck et al.ll2005l : lYan et al.ll2005l : iLutz et al. 
20051 ). and typically sample only the most extreme ULIRGs (> 10^^ Lq). 



MIR spectra of starbursts are dominated by emission from polycyclic aromatic hydro- 
carbons (PAHs) at 6.2, 7.7, 8.6, 11.3, 12.7, an d 17 /xm and broad (~ 2 /xm) silicate absorption 
at 9.7 and 18 fim (iGenzel fc CesarskyI |2000| . and the references therein). In local sources. 
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the presence of PAH features, combined with the shape of the continuum, can be used to 
determine the nature of the emi ssion (AGN vs . star formation) that don i inates the bolomet- 



ric luminosity of a source (e.g. iLaurent et aLll2000l : iPeeters et aLll2004l : lArmus et aLl 12006 



Spoon et aLll2007l ). The rehabihty of these same diagnostics at high redshift requires further 
demonstration. 

LIRGs at 2 ~ 1 and ULIRGs at 2; ~ 2 are hkely to have flux densities on the o rder 
of 0.1 mJy at ~ 20 /im, based on the bolometric corrections of IChary fc ElbazI (120011 ). In 
principle, such sou rces should be observable with the Spitzer Infrared Spectrograph (IRS; 
Houck et al.l 120041 ) in several hours, assuming that the noise characteristics do not change 
at faint fluxes and the signal-to-noise ratio scales with the square root of observing time. 
In order to push the limits of MIR spectroscopy with Spitzer, we targeted two faint (0.15 
mJy) galaxies at 2; > 1, using 32 hours of Director's discretionary time. In this paper, we 
present the results of those observations, demonstrating for the flrst time the capability of 
the IRS to detect continuum emission at such faint flux levels. Throughout, we assume a 
A-dominated flat universe, with Hq = 71 km s~^ Mpc~^, = 0.73, and Qm = 0.27. 



2. Target Selection 



We selected targets from the southern GOODS fleld. GOODS, by virtue of being 
the target of the deepest Spitzer, Hubble, Chandra, and ground-based imaging and spec- 
troscopy provides the richest, most homogeneous dataset yet compiled to understand the 
formation and evolution of galaxies and AGN. While the Northern fleld offers certain advan- 
tages, the Southern fleld was better positioned for the available Spitzer scheduling window. 
Amon gst the unique compo nents of the Southern fleld are: the Hubble Ultra-Deep Field 
(UDF; iBeckwith et al.ll2006l). which pro vides high resoluti on ACS imaging; extensive spec- 



troscopic redshifts (ILe Fevre et al. 



shallow 16 /im photometry (ITeplitz et al.ll2005l ): and ultradeep (1 Ms) Ch andra 2-10 keV 



20041): ultrad eep 24 fj,m (IChary et al.ll2006. in prep.l ) and 



X-ray observations which provide the most complete census of AGN activity (iGiacconi et al. 
2OO2I ). 



The goal of the observations was to examine relatively faint IR-luminous galaxies at 
high redshift. We estimate the IR luminosity from the obser ved 24 /im flux. The MIR 
luminosity of local galaxies in the IRAS bright galaxy sample ( ISoifer et al.l 119871 ) has been 
found to correla te strongly with their far-infrared luminosity, which is dominated by large, 
cool dust grains (jchary &: Elbaz 2001 ). This correlation has been applied to develop a library 
of model templates of the mid- and far-infrared SEDs of galaxies. The library consists of 
template SEDs across a range of luminosities, which can be redshifted to predict the MIR 
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flux of a source with given luminosity at a redshift of interest. For each source in UDF, 
we select the template for which the library predicts the closest 24 /im flux density at the 
appropriate redshift to apply a bolometric correcti on; we do not use the shor ter wavelength 
IRAC measurements. The corrections based on the lChary &: Elbaa (120011 ) and lDale &: Helou 
(I2OO2I ) template are used to derive an infrared luminosity (Lir= 8 — 1000 fim). The use of 
/24 to predict Ljr using these templat es is suppo r ted by the observation of sub-millimeter 
galaxies (SMGs) in GOODS-North bv lPope et aP J2006h . 



We chose two sources in the UDF, one with a Chandra detection and one without. 
We chose each source to have a ratio (or upper limi t) of X-ray to infrared luminosity of 
Lx/LiR< 1.5 X 10~^, consistent with many SMGs ([Alexander et al.ll2005l ). We avoided 
sources that are only marginally resolved in HST imaging, which might indicate strong 
AGN contributions, and selected instead from those with fully resolved morphologies. A 
final criterion was the spectroscopically-confirmed redshift, which we selected to place the 
strongest PAH features in the center of one of the IRS passbands. By selecting one source 
at 2 ~ 1 and one at 2; ~ 2.5, we were able to obtain spectra covering the 6.2, 7.7 and 8.6 /im 
PAH features in just one or two slits. 

To meet these criteria, we selected a. z = 2.69 C/ianc?ra-detected source with 24 /im flux 
density of ~ 0.13 mJy, hereafter Source 1, and a z=1.09 object with a 16 fim flux density 
of ~ 0.15 mJy and no X-ray detection, hereafter Source 2. The infra red luminosity o f these 
sources is approximately 10^^ and 3 x 10^^, respectively, based on the lChary &: Elbaa (1200 ll ) 
templates. The Chandra detection of Source 1 is only in the soft-band (0.5-2 keV), which 
does not immediately rule out star-formation as the origin of the X-ray emission. Given the 
experimental nature of the observations, we were given more than usual latitude to arrange 
the precise scheduling of the observations. Thus, we were able to choose observation dates 
which were favorable for including additional sources in the slit. Figure [T] shows the position 
of the slit overlayed on the 16 fim map of the field. Two of the secondary objects in the 
slits are of great interest, and are indicated in the figure: (1) an X-ray source detected by 
Chandra in the soft-band only, but lacking an optical redshift, referred to hereafter as Source 
2-x; and (2) a source lac king an optical red shift but selected to be at high redshift by the 
52;/^ selection technique (jPaddi et al.l 120041 ). referred to hereafter as Source 1-BzK. Figure 
12] shows the ACS images of each source. Table [T] gives the photometry for the four sources 
of interest. Both LL slits also include a somewhat brighter object (/le ~ 0.25 mJy) which 
appears to be a single source in the 16 fim map, but IRAC imaging reveals it to be a blend 
of three sources. The slits were oriented to include this source in case the sensitivity proved 
worsen than predicted and the primary targets were not detected. Given its blended nature, 
it is not discussed further in the present paper. 
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The 16 /iin map shows that Sources 1 and 1-BzK have close, but clearly separated, 
companions of comparable brightness. They lie ~ 8" from the targets, which is less than 
two pixels in LL-1 (5.1" pixels). Although we do not know that the companions lie at 
different redshifts than the targets, we have no evidence that they do; we expect that they 
will not affect the detection of emission features, though they will make the measurement 
of equivalent widths less certain. The companion sources are of comparable brightness to 
the targets, so, if half of the companion light were to fall into the extraction aperture, the 
continuum would be overestimated by about 50%. 

3. Observations and Data Reduction 

We chose integration times that would likely achieve signal to noise ratios of 2 — 4 in the 
continuum, so that the results could be evaluated even if no emission or absorption features 
were detected. We obtained the observations using a different observing mode {Staring or 
Mapping) for each object, in order to compare the effectiveness of the modes. 

For the z = 1.09 target. Source 2, we obtained a spectrum using the second order of the 
Long-Low module (LL-2, 14 — 21 /im). The spectrum was taken in Staring Mode, in which 
the target is placed at two "nod" positions within the slit. The observation was broken into 
two separate Astronomical Observing Requests (AORs), each of which contained 4.5 hours 
of on-source integration divided between the two nod positions. At each nod position we 
took 70 spectra of 120 sec ramp duration. The total on-source integration time was 9.4 
hours, and the total AOR time with overhead was 11.9 hours. We also observed this target 
in the short- wavelength low-resolution slit, in the first order (SL-1), in a single AOR for 
three hours on source in Staring Mode. 

For the higher redshift target. Source 1, we used the first order of the Long-Low module 
(LL-1, 21 — 37 iim). These observations were obtained in Mapping Mode, wherein the target 
is placed at multiple positions along the slit; we selected six map positions. The observations 
were broken into three AORs, each with 20 spectra of 120 sec ramp duration at each of the 
six map positions. The total on-source integration time of 12 hours, and the total AOR time 
with overhead was 15.4 hours. The observations are summarized in Table [2l 

Observations were obtained on 06 August 2005 (SL-1), 10 August 2005 (LL-2), and 12 
September 2005 (LL-1). The LL-2 AORs were executed consecutively. Two of the LL-1 
AORs were consecutive, but there was a gap of more than 12 hours between the second 
and third observation. We observe no degradation in signal to noise resulting from the gap. 
All observations were scheduled immediately after the "skydark" calibrations in order to 
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maximize the sensitivity by reducing the change of latent images from a preceding bright 
target. 

To guard against the unhkely event of large deviations in telescope pointing, we obtained 
high accuracy peak-ups in the blue filter, using a star selected from the 2MASS catalog. 
During observation planning, we examined the MIPS 24 /im image of the star prior to 
submitting the AORs to ensure that its flux was accurately predicted and that it had no 
bright neighbors. We inspected the PU data, to check that the centroid was reasonable and 
that the pointing was as expected. No problems were identified. 

3.1. Data Reduction 

The Basic Calibrated Data (BCD) were produced by the S13 pipeline at the Spitzer 
Science Center (SSC), which includes ramp fitting, dark sky subtraction, droop correction, 
linearity correction, fiat fielding, and wavelength calibration (see the IRS Pipeline Hand- 
boolll]). Further processing of the two dimensional dispersed frames is required before spec- 
tral extraction. In addition to the pipeline processing, we performed the following reductions: 
latent charge removal, rogue pixel interpolation, and residual sky subtraction. 

A small fraction (1-2%) of charge on the detector persists between frames despite the 
resetting of the detector that occurs prior to every integration. This latent charge decays 
slowly over time and is removed by the annealing process. In the case of very faint sources, 
the source of latent charge is the zodiacal background. Over the course of a six hour AOR, 
this charge can build up to a significant level. We removed the latent signal by fitting the 
slope of the background with time. We measured the median background independently for 
each row of each frame, fitted the slope of these median values, and then subtracted this fit 
row-by-row from each frame. The background in BCD images is the residual after "skydark" 
subtraction in 3D. In the deep observations, the residual sky level was ~ 50 electrons per 
second, and the latent charge built up to ~ 5 electrons per second by the end of the AOR. 

Unstable, or "rogue" , pixels are those which are usable in some AORs and not in others, 
depending on the recent history of the detector. A mask of known rogue pixels is provided 
by the SSC, and we identified further suspect pixels from the data themselves. We used the 
IRSCLEAN program provided by SSC to find rogue pixels in the 2D data, using its default 
settings. We further searched for pixels with abnormally high variance (> lOo") with time. 
All pixels identified as rogue were interpolated over using IRSCLEAN. 



http: / /ssc. spitzer. caltech.edu/irs/dh/ 
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Once the BCD data were "cleaned" , we created residual sky images from the data. In 
staring mode, wc used the other nod position, and in mapping mode we used the other five 
map positions. In both cases, we used a resistant mean (removing outliers before calculating 
the mean) with time to calculate the sky value in each pixel. 

The individual reduced frames were coadded to produce final 2D spectra at each nod or 
map position. One dimensional spectra were extracted using the SPICE software provided 
by the SSC. To minimize the noise contribution from the background, we used narrower 
extraction windows than the default. We used a window which expanded with wavelength 
and had a width of ~ 2 pixels at the blue end of each order. The "slit losses" introduced by 
the narrow window were estimated by extracting the spectrum of the standard star HR 7341 
(taken in the same campaign) using both the standard and the new width. The correction 
factor varied with wavelength by ~ 10% and had a mean value of 1.4. The spectra of 
HR 7341 were taken at the standard Staring Mode nod positions, so the Mapping Mode 
spectra may not be perfectly corrected (at the < 5% level), due to fiat field variation or pixel 
undersampling. 



3.2. Noise Measurement 

The signal to noise ratio (SNR) in the IRS spectra was measured using the variance 
of each pixel from frame to frame. We measured the variance in the sky frame for each 
pixel, scaled by e.g. a/2 for the nodded frames, and created two dimensional noise frames 
for input into SPICE. The uncertainties were then propagated by the extraction process. 
These uncertainties were checked against the standard deviation in the continuum of the ID 
spectra at, for example, ~ 24 /xm in LL-1. 

The SNR at 16 in LL-2 and 24 fim in LL-1 is generally consistent with the predictions 
of the SSC website (SPEC-PET). The predictions do not include noise introduced by the 
pipeline processing or spectral extraction. They also assume the spectra are smoothed to a 
resolution of R = 50, which we have not done. A 12 hour integration in LL-1 is expected to 
yield a 3cr uncertainty of 0.12 mJy at 24 /im with spectral resolution smoothed to R = 50. 
Our measured 3a uncertainty is ~ 0.13 mJy. SPEC-PET suggests a 9 hour integration 
in LL-2 will yield a 3a uncertainty of 0.08 mJy at 16 /im and R — 50. Our measured 
uncertainty is ~ 0.12 mJy. While both spectra achieve sensitivity close to the estimates, the 
Mapping Mode observation gives noticeably better results. Many factors may contribute to 
the difference between the modes, including better sky estimation in Mapping Mode, as well 
as an averaging over low level variation in the fiat field. 
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The SNR varies strongly with wavelength, with noise rising sharply at the red end of 
each order. The shape of the error array measured in our data agrees well with the predictions 
of SPEC-PET, independent of the overall level. 

4. Results 

Figure [3] shows the 2D spectra of the sources. Both primary sources and secondary 
targets discussed above (Source 1-BzK and Source 2-x) are detected. Figures HHT] show the 
extracted spectra of the primary and secondary targets. The continuum level is in good 
agreement with the estimate from the broad-band photometry at 16 and 24 fim. Broad 
emission features, identified as polycyclic aromatic hydrocarbon (PAH) emission, are de- 
tected at rest wavelengths of 6.2, 7.7, and 8.6 /im, along with the continuum at the blue end 
of LL-2 and LL-1 (except in Source 2-x). 

The PAH features in Source 1-BzK confirm that it is at high redshift, indicating z = 
2.55 lb 0.03. The uncertainty here is based upon the narrower 6.2 /im PAH emission feature, 
and excludes the broader 7.7 fim feature. This redshift is consistent with the expectation 
from the optical photometry. The other secondary source. Source 2-x, shows a single definite 
emission line at an observed wavelength of Xobs ~ 16 /im, and possibly broad enhanced 
emission at Xobs ~ 19 /xm. If the isolated line were identified with the 6.2 fim PAH feature 
and the broad emission with the 7.7 /im PAH feature, the object would lie at 2; ~ 1.55,. 
Identification of the isolated emission line as the 7.7 /xm PAH feature (indicated 2; ~ 1.1) is 
less likely, as that would center the broad emission at rest-frame 9.5 fim, which might only 
be explained by strong silicate emission, which is unusual in a source with PAH emission. 

Many infrared luminous sources show deep absorption troughs at rest-frame ~ 10 /xm 
from silicates. This rest-frame wavelength is at the red end of the order for all target spectra 
and thus the presence of the absorption is difficult to establish. 

We present estimates of the equivalent width (EQW) of the PAH emission lines in each 
object in Table [31 The continuum was measured on both sides of the line, when detected, 
and assumed to have negligible slope across the line. The continuum in all four objects is 
similar, showing shallow or no slope across the modest restframe wavelength range covered 
by each order. 
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Discussion 



Sources 1, 2, and 1-BzK all show strong PAH emission. Source 2-x has emission that 
may arise from PAHs as well. Figur e [8] com pares the spectra of each source with the starburst 
template spectrum of iBrandl et al.l (120061 ). The PAH features in all four sources are broadly 
similar to the template. Local starbursts and star-forming ULIRGs typically have rest- 
frame equivalent widths between ~ 0.5 and ~ 0.8 /im for both the 6.2 a nd 7.7 /im features , 



and equivalent wi dth ratios between these features of around unity (IBrandl et al.l 12006 



Armus et al.ll2006l ). Our UDF sources show slightly lower rest-frame EQW(6.2), and the 
ratio of EQW(7.7)/EQW(6.2) appears larger by a factor of about two in the z > 2 sources 
and a factor of about 1.5 for the z = 1.09 source. This difference is largely attributable to 
the low continuum near the 7.7 /xm feature in the high redshift objects. The continuum is 
detected only at 3-4 a confidence, so this result is hard to quantify. The presence of close 
companions to the z > 2 sources also contributes to the lower EQW (see section [2]). 

At low redshift, PAH emission is usually associated with star formation. Sources with 
strong AGN generally have very weak or absent PAH emission. The strength of the PAH 
emission can thus be used to estimate the relative contribution of star formation to the 



bolometric luminosity (e.g. iLaurent et al.l l2000l : iPeeters et al.l 120041 ). so these spectra are 
suggestive of starburst sources. It is not surprising that Source 1-BzK or Source 2, a z = 1.09 
LIRG, show prominent PAH emission. However, one might have expected Source 1, which 
is detected in X-rays, to have weaker emission. 

The wavelength coverage of the present observations is quite limited. We can extend 
the SED by considering additional observations of the GOODS fields: Spitzer IRAC data at 
3.6, 4.5, 5.8, and 8 /im; the 16 /im IRS pho tometry (Teplitz et al. 2007, in prep.); and the 
HST/ACS and ground-based NIR imaging (iGiavalisco et al.l 120041 ). These extended SEDs 
are shown in Figures [H The figure inset shows the wavelength region around the ~ 1.6 
/im opacity minimum in the stellar atmospheres of AGB stars. This photospheric emission 
"bump" is commonly observed in galaxies with large evolved stellar popul ations and is often 
used to identify potential massive starbursts (e.g. iLonsdale et al.ll2005l ). as it is sensitive 
to a lack of hot dust emission, which would drown out the "bump." However, the feature 
may also be prominent in relatively modest star-forming galaxies if they lack significant dust 
emission. 

The extended SED shows that Source 2, at 2; = 1.09, appears indeed to be a typical, 
stellar-dominated LIRG, as predicted by its other properties. A clear emission peak is seen 
in the IRAC photometry at 1.6 /im in the rest-frame. 



Source 1 ai z = 2.69, on the other hand, does not show this "bump." In fact, its 
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SED appears relatively featureless at all wavelengths shortward of 6.2 /im, extending into 
the NIR and optical. This featureless SED is consistent with the expectation for an X-ray 
source at high redshift. Source 1, therefore, is likely to be a composite object containing an 
AGN which dominates the emission in the near-IR, swamping the 1.6 fim bump, together 
with significant ongoing star formation which is responsible for the PAH emission in the 
rest-frame MIR. We also find that source 1-BzK has a MIR SED similar to that of Source 
1, lacking a prominent 1.6 fim bump. 

If Source 1 is, indeed, a composite object, then our estimate of its luminosity, ~ 10^^ L0 
based on a starburst template, may be too high. We can refine this estimate by scaling the 
luminosity of a local composite source. UGC 5101 is a local ULIRG that shows significant 
PAH emission, but also has optical line ratios consistent with an AGN which contributes 
less than 10% of the bolometric l uminosity, but dee p silicate absorption leaves open the 
possibility of a larger contribution (lArmus et al.ll2004l ). Normalizing the SED of UGC 5101 
at 24 yum to Source 1, we estimate the bolometric luminosity of Source 1 to be 4 x 10 Lq, 
a factor of about 2.5 less than using a pure starburst template. If we attribute all of the 
flux in Sourc e 1 at restframe 5 /xm to an AGN, and scale a representative AGN template 



(05189-2524: lArmus et al.ll2006l ) to that flux, we find a contribution of 1.6 x 10^^ Lq from 
the AGN to Ljr. 

Using the UGC 5101 template suggests ratio of Lx/-Z^ir~ 0.004 in Source 1, consis- 
tent with the ratio ob served for radio-detected SMGs with significant AGN contribution 
( lAlexander et al.ll2005l ). The C/iandra non-detection of Source 1-BzK and Source 2 sets lim- 
its on their Lx/Lyr ratios of about 0.001, which does not rule out the presence of a weak AGN 
in those objects. Source 2-x has a ratio of ~ 0.007 (similar to that of an AGN-dominated 
SMG) if it is at the adopted redshift. 



5.1. The 1.6 fim Bump vs. PAH emission 

Source 1 likely includes a significant AGN contribution, combined with strong PAH 
features but lacks a 1.6 /xm bump. To quantify the sensitivity of the bump to the strength 
of an AGN contribution, we construct a simple model of composite objects by mixing the 
SEDs of a local starburst (NGC 7714) and a pure QSO (PC 0804+761), both of which are 
first normalized to have unit 1-1000 fim luminosity. We find that even a weak (10% of 
the mixture) AGN may produce enough hot dust emission in the NIR to hide the 1.6 /im 
bump, while leaving the PAH features easily detectable. This hot dust continuum does, 
however, noticeably reduce the EQW of the PAH emission. We plot a 10% AGN composite 
in comparison to Source 1 in Figure [SI 



We can compare the PAH emission with the stellar bump in a variety of other objects. 
We can quantify the strength of the 1.6 fim bump using the ratio of the rest-frame near-IR 
J/K to K/H slopes: 



B 



[W)/m)] 

[UK)/UH)] 



(1) 



This parameter is large when the slopes show an inflection, as expected for a strong 1.6 
/im bump, but it is small when the slope of the SED is constant. In Figure [10] we show 
the strength of the B parameter versus the equivalent width of the 6.2 /xm PAH emission 
feature for the UDF sources compared to the local ULIRGs and other high redshift targets. 
The value of B was obtained by interpolating the observed SED to get the rest-frame JHK 
photometry. In the figure, we indicate the average strength of the bump in a sample selected 
on the basis of that f eature, by plotting the average value of B for objects observed by 
Weedman et al.l (120061 ) . who obtained IRS spectra of sources in SWIRE fields with strong 
1.6 /im bumps. Note that the plotted value is not a selection limit, merely an indicative 



value. The local sources used for comparison are local U LIRGS observed bv lArmus et al. 
(120061 ). The near-I R photometry is tak en from 2MASS (ISkrutskie et al.ll2006l ). using the 



Large Galaxy Atlas (jjarrett et al.ll2003l ) where available and the Extended Objects catalog 
otherwise. We also plot the ratio that would be expected given the NGC 7714 plus QSO 
model for composite objects. We find that there is a general trend of high PAH EQW for 
objects with large 1.6 yum bumps. While the uncertainties are large, the trend is about what 
is expected from the composite spectrum. 

A number of factors complicate Figure [TOl First, the model assumes that there is little 
extinction at 1.6 /xm. Strong dust reddening will reduce the value of B by lowering the J- 
band flux relative to the ii'-band. The effect of additional reddening is shown in the figure. 
Secondly, the strength of the bump depends on the fraction of old stars in the source while 
the PAH emission depends on the number of young stars. In a sense, the average age of 
the population and the amount of dust determines an object's location in Figure [TOl Thus 
a source with a much older underlying population than NGC 7714 would have a higher 
ratio of B to EQW(6.2) than the composite model. Conversely, a very young object lacking 
significant old stars could have a much lower ratio. Finally, we have implicitly assumed that 
the geometry of all sources is similar. Both the Spitzer photometry and spectroscopy of 
high redshift sources capture the light of the entire source, because high redshift sources are 
unresolved and fall within the entrance slit of the IRS. For lo cal ULIRGs, on the other hand, 
the IRS entrance slit usu ally includes most of the MIR flux ( Soifer et al. 200ol: Armus et al. 



20061 : ISurace et al.ll2007l ). but is much smaller than the near- IR extent JScoville et al.lboooh . 
Therefore, for the local ULIRGs in Fig. [10] we have used the total NIR light instead of an 
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aperture matched to the size of the IRS sht, simply to make comparison to the high-z sources 
less problematic. 



5.2. Comparison with Other Surveys 



The set of high redshift IRS spectra considered here is small, but we examine if it is 
typical of other objects targeted in this flux range. 

First, we compare the properties of the faint, high-z starbursts with more luminous high- 
z objects previously observed with the IRS. The largest samples were selected either using 
optical-to-MIR co l or properties or by s election at a longer wavelength, e.g. the sub-millimeter 
Jlutz et allboosl l | Houck et al.l (120051 ) selected objects based o n very faint { Iar > 25) optical 
magnitudes and detectable (> 0.8 mJy) 24 /xm flux densities. lYan et al.l (120051 ) used a less 
severe R — [24] selection function and also included the requirement of blue 8 to 24 fim color. 

The sources in the present sample are brighter in the optical than the Houck et al. 
sources, being detected at Iab ~ 24.5 while the latter were undetected at the same brightness. 
The UDF sources are much bluer, as the Houck et al. sources are 5-10 times brighter at 24 
/im. This difference is not surprising given the rarity of the Houck et al. sources (33 in 8 
square degrees). The z > 2 sources in the present sample meet the Yan et al. selection 
criteria, but the lower redshift sources are bluer in — 24 color. The Houck et al. and 
Yan et al. samples find substantial numbers of objects with spectra showing little PAH 
emission, likely indicating AGN-dominated sources. These objects are brighter than those 
in the current data by a factor of ~ 5 — 7, suggesting bolometric luminosities in excess of 
10^^ L0 if starburst-dominated templates are used; if the sources were AGN-dominated, a 
smaller bolometric correction would be appropriate and they might be closer in luminosity 
to the UDF sources. Two sources in the Yan et al. sample have measured PAH emission; the 
EQW(6.2) in one object is consistent with the UDF sources while the other has significantly 
weaker PAH emission. 

The GOODS-South region has no t been obs e rved at submm wavelengths, so a direct 
comparison with SMGs is not possible. iLutz et al.l (120051 ) obtained IRS spectra of two SMGs 
at z ~ 2.8, one of which shows strong PAH emission while the other is dominated by a red 
AGN continuum over the same rest-frame wavelength region covered by the UDF spectra. 
Lu tz et al. place the AGN - dominated SMG in the same region of the Lx / -Z^ir plane identified 
by lAlexander et al.l (120051 ). iMenendez-Delmestre et al.l (12006! ) have also observed SMGs with 
the IRS. They find the SMG spectra to be consistent with that of a composite of a strongly 
star- forming galaxy and weak AGN. The EQW of the 7.7 /xm PAH feature in their spectra 
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is consistent with that of the z > 2 UDF sources. 

As discussed above, Source 1 at 2; = 2.69 provides a striking example of an object that 
would be excluded by a "bump" survey but nonetheless may have star formation contribut- 
ing the majority of its bolometric luminosity. Sources of this kind may be an important 
contribution to the cosmic star formation density. However, based on current data, it is dif- 
ficult to estimate how many such sources will be missed in IRAC surveys. This uncertainty 
demonstrates the need for wider use of deep MIR spectroscopy of high redshift sources. 



6. Summary 

We have obtained the longest exposure-time spectra to date using the the Spitzer IRS, 
demonstrating the power of the IRS to trace the properties of LIRGs/ULIRGs at z > 1. 
The spectra of four high redshift sources in the Hubble UDF show prominent PAH emission, 
comparable to local starbursts. The combination of IRAC photometry, IRS spectroscopy and 
X-ray data on these sources place constraints on the physical properties which dominate their 
luminosities. These sources follow the rough correlation in ULIRGs between the strength of 
PAH emission (arising from young stars) and the 1.6 /xm bump (arising from older stars). 



This work is based in part on observations made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. Support for this work was provided by NASA through an award issued by 
JPL/Caltech. 
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Table 1. Photometry 



Obj. 


R.A. 


Dec. 


Iab 


/8 


/l6 


/24 




(J2000) 


(J2000) 


(mag.) 


(mJy) 


(mJy) 


(mJy) 


1 


3:32:44.00 


-27:46:35.0 


24.47 


0.0201 ±0.0004 


0.05 ±0.025 


0.126 ±0.010 


1-BzK 


3:32:38.52 


-27:46:33.5 


24.65 


0.010 ±0.005 


0.10 ±0.025 


0.200 ±0.020 


2 


3:32:34.85 


-27:46:40.4 


23.30 


0.0153 ±0.0004 


0.15 ±0.025 


0.133 ±0.015 


2-x 


3:32:39.72 


-27:46:11.3 


24.96 


0.0467 ± 0.0004 


0.16 ±0.025 


0.225 ±0.020 
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Table 2. Observations 



Obj. 


z 


Sub-slit 


Mode 


Ramp 
Time (s) 


Cycles 
per AOR 


AORs 


1 


2.69 


LL-1 


6x1 map 


120 


20 X 6 pos. 


3 


2 


1.09 


SL-1 


Staring 


240 


25 X 2 nods 


1 






LL-2 


Staring 


120 


70 X 2 nods 


2 



Table 3. Emission Features 



Obj. 


z 


log(LiR/Lo)'^ 


Flux(6.2)^ 


(EQWrest)^ 


Flux(7.7)^ 


(EQWrest)^ 


Flux(8.6)^ 


(EQWrest)^ 


1 


2.69 


12.6 


5.7 ±0.7 


(0.24 ±0.06) 


23.3 ±2.0 


(0.98 ±0.23) 


5.7 ±0.7 


(0.39 ±0.09) 


1-BzK 


2.55 


12.6 


6.4 ±0.5 


(0.24 ±0.09) 


20.1 ±2.1 


(0.79 ±0.29) 






2 


1.09 


11.5 


15.3 ± 1.0 


(0.6 ±0.2)^ 


29.9 ± 1.6 


(0.90 ±0.22) 


5.8 ± 1.6 


(0.22 ±0.07) 


2-x 


~ 1.55^ 


12.3 















^The infrared luminosity inferred from the use of local templates as discussed in Section [51 As only one local 
template is used per o bject, the uncertainty is not given. Typical uncertainties using a suite of templates is < 40% 
Jcharv fc ElbazlboOlh . 

^lO""*^^ ergs cm~^ s~"^ S 

■^EW given in /im. In most cases, the uncertainty is dominated by the uncertainty in the continuum. 

*^The continuum level is only marginally detected in the SL spectrum, so this measurement is uncertain 

''The redshift estimate is based on the identification of an isolated emission line as the 6.2 /im PAH line. The 
continuum is not well detected, so EQW values are not given. 
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Fig. 1. — The orientation of the LL-1 (red) and LL-2 (yellow) slits overlayed on the 16 //m 
image of the UDF. The image size is 4' x 2.5'. The primary sources are indicated by their 
redshift, as are other sources in the slits. 
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Fig. 2. — Three-color ACS images of the four sources of interest. The images are the 
combination of the Viz filters in the GOODS data. Images were obtained from the MAST 
high level science products. The images are 5 arcseconds on each side; north is up in the 
images, and east is to the left. 
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Fig. 3. — The 2D dispersed frame image of the IRS spectra of faint sources. The top panel 
shows the LL-1 spectrum of Source 1 and Source 1-BzK, in one of the six map positions. 
The bottom panel shows the LL-2 spectrum of Source 2 and Source 2-x in one of the two 
nod positions. Many of the pixels just above Source 2-x have been masked. In each case, 
only the illuminated portion of the relevant order is shown. The approximate wavelength 
scale is indicated. Positive signal is black on white. 
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Fig. 4. — The spectrum of Source 1 aX z — 2.69 in LL-1. The la error arrays are shown 
(dashed hnes). 
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Fig. 5. — The spectrum of Source 2 dX z = 1.09 in SL-1 (thin hne) and LL-2 (thick hne). 
The 1(7 error arrays are shown (dashed hues). 
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Fig. 6. — The spectrum of Source 1-BzK a,t z — 2.55 in LL-1. The la error arrays are shown 
(dashed hues). 
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Fig. 7. — The spectrum of Source 2-x in LL-2, likely at z ~ 1.55. The la error arrays are 
shown (dashed hues). 
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Fi g. 8. — The spectra of each of the four UDF sources compared to the starburst template 



of iBrandl et al.l (120061 ) (smooth sohd hne). The UDF source spectra have been smoothed by 
2 wavelength elements. The template spectrum has been normalized to match the peak of 
the PAH emission features. The spectrum of Source 1 is also compared to the model of a 
starburst- AGN composite (dashed line; see Section [5?T1) . with a 10% AGN contribution. 
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Fig. 9. — The spectral energy distributions (SEDs) of the four UDF sources, using the the 
IRS spectra together with the HST (triangles), ground-based NIR (squares), and Spitzer 
(circles) photometry. Inset plots show the wavelength region around 1.6 /im (dotted line) on 
a linear scale. 



- 29 - 



0.8 



X 0.6 

ffi 0.4 

< 

J 0.2 

a 



0.0 



SWIRE 

average 



,'x=0 



10%/' 



50% 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

[f/fK]/[fK/fH] 



Fig. 10. — The strength of the 1.6 /iin bump, parameterized as the ratio of NIR slopes, vs. 
the rest-frame EQW of the 6.2 //.m PAH fe ature (in micron s). We compare UDF sources 
(circles; this work ), FLS sources (triangles: lYan et al.ll2005l ). and local ULIRGS (squares; 
Armus et al1l2006h . For local ULIRGs, we have used the total NIR light to measure the 



bump, rather than an aperture matched to the size o f the IRS slit. Th e dotte d line indicates 
the average strength of the bump in the sample of IWeedman et al.l (120061 ) , selected from 
the SWIRE survey. The dashed curves indicate the expected locus for the mixture of a 
starburst and a quasar (see Section 15. ip . The r = curve has no additional reddening 
beyond that in the starburst spectrum, while the r = 1 curve includes an addition reddening 
of r(9.7yum) = 1. Labeled tick marks indicate the 10% and 50% AGN contribution to the 
mixture of normalized 1-1000 /xm luminosities. 



